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Abstract: Non-Abelian gauge theories may have continuum limits in more than four 
dimensions, supported by non-trivial ultra-violet fixed points. Moreover, such theories 
can be expected to be accessible to Wilson's epsilon expansion. We investigate this series 
for SU (N) Yang-Mills, in particular for the fixed point coupling and exponent up to 
four loops. From the model-building point of view, such theories would be effectively per- 
turbatively renormalizable in the normal way. A particularly attractive possibility is the 
construction of renormalizable extra-dimensional models of the weak interactions, which 
have the potential to address the full hierarchy problem. The simplest such gauge-Higgs 
unification model is however ruled out by a combination of theoretical and phenomenolog- 
ical constraints. 
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1. Introduction 

There has been considerable interest in the last few years in constructing field theories of 
particle physics in more than four dimensions. On the one hand, these ideas are inspired by 
similar features in string theory, on the other hand, elegant new approaches to the persistent 
mysteries of theoretical particle physics are now possible [1, 2]. These aim variously to solve 
the famous hierarchy problem (i.e. the fact that the Higgs mass is quadratically sensitive to 
higher scales), the little hierarchy problem (the fact that precision measurements seem to 
favour a scenario with a low Higgs mass ~ 115 GeV and no new physics until ~ 10 TeV [3]), 
the pattern of fermion masses and mixing angles, and problems with grand unified theories 
such as doublet-triplet splitting, proton lifetime and mass relations. The new approaches 
to solving these issues include theories with so-called universal extra dimensions which e.g. 
result in calculable SUSY breaking and Higgs masses [4] , gauge-Higgs unification in which 
the Higgs field is a component of a higher dimensional gauge field whose mass is thus 
protected by gauge invariance [5-13], Higgsless models where the symmetry is broken by 
the boundary conditions [14], grand unified models in which colour triplet Higgs are given 
a Kaluza-Klein mass whilst the doublet remains massless at the compactification scale [15] 
(see also [16]), and approaches to the flavour problem using the freedom to place fermions 
at different points (branes) in the extra dimensions and to couple them non-locally through 
other fields or Wilson lines [9, 13, 17-20]. 

However all of these new approaches suffer from a severe drawback: quantum field 
theory as conventionally envisaged is not renormalizable in more than four space-time di- 
mensions. Thus these theories must be viewed as effective field theories, only applicable 
over a limited energy range. For energies much lower than the defining scale E in the 
theory (e.g. E ~ 1/R where R is the compactification scale) experimental constraints 
require that it reduce to the Standard Model plus very small corrections. On the other 
hand there is a maximum energy A, usually taken to be the energy where some coupling 
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becomes strong and estimated by so-called naive dimensional analysis (NDA) [21], above 
which something other than quantum field theory must take over (for example some conjec- 
tural string theory). This maximum energy depends on the couplings, with larger A/E for 
smaller couplings, but typically A ~ 10E to 100-E is possible. With the ultraviolet comple- 
tion unknown, we can only bound the symmetry-preserving higher dimensional vertices in 
the effective theory, under the assumption that they naturally have order one coefficients at 
the cutoff scale. This sets an irreducible limit on the expected predictivity of these models 
controlled by E/A, again dependent on the details, but typically authors aim for ~ 1%. 

This limitation is less of a problem for theories that attempt to describe the energy 
range between the GUT scale 10 15 GeV and the Planck mass 10 19 GeV, and for which the 
irreducible uncertainty of ~ 1% comes on top of already (better be) small corrections to 
Standard Model quantities. 

It is a much more drastic issue for higher dimensional theories that directly address 
the hierarchy problem. These models tackle the problem by introducing new physics at 
E ~ 1 TeV in such a way that the Higgs mass becomes insensitive to higher scales. Here 
attempts to bring such a model into agreement with the impressive precision LEPI/II 
data can in turn force it into regions of parameter space where the need for large ratios 
in couplings then result by NDA in significant reductions in A, cutting off its domain of 
applicability and thus severely limiting its predictability (e.g. see [13,22]). Since an as yet 
unknown ultraviolet completion at some intermediate scale A is required, these models can 
actually only address the little hierarchy problem, because any insensitivity of the Higgs 
mass to scales higher than A has to be explained in the unknown theory. It is in any case 
a rather drastic step to address the little hierarchy problem by abandoning quantum field 
theory altogether at say ~ 100 TeV for a largely conjectural ultraviolet non-field theory. 

It might be possible however, to rest certain restricted classes of such models on a 
more secure foundation. In the 1970s it seems to have been commonly understood that 
non-Abelian Yang-Mills theory may exist as a non-perturbative quantum field theory, a.k.a. 
continuum limit, in more than four space-time dimensions [23]. With few exceptions since 
then [24-31] this possibility seems to have been ignored or forgotten. Independently in 
string theory, probably related supersymmetric versions have been discovered [32]. 1 

We briefly review the little evidence that has collected, for and against such theories, 
in sec. |2| We also investigate the series for the fixed point coupling d*, and the critical 
exponent u, using the e expansion up to four loops. Together, the evidence suggests that 
these fixed points exist in D = 5 dimensions, possibly in D = 6 dimensions and probably 
not for D > 7 dimensions. 

At any rate such theories formally do exist within an expansion in e of a 4 + e di- 
mensional theory. As we review in sec. §, this means that quantum corrections can be 
computed perturbatively and renormalized in four dimensions. From the model build- 
ing point of view, the theory is perturbatively renormalizable; one only needs to bear in 
mind that at high energies the appropriate dimensionless ratio d, involving the Yang-Mills 
coupling a, runs to a computable finite value rather than zero. 

1 A Kaluza-Klein treatment appears in ref. [33]. 
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For the continuum limit to exist within the e expansion when other fields are included, 
one requires that in four dimensions, all couplings are relevant (effectively masses and 
couplings of positive mass dimension) or marginally relevant (equivalently asymptotically 
free). We outline the restrictions this imposes in sec. |3[ We then extend this understanding 
of renormalizability to models with branes, and discuss the constraints that arise, concen- 
trating on the bosonic sector, as well as making some further comments on the applicability 
of this idea to Grand Unified Theories. 

A particularly attractive application is to extra dimensional models of the weak inter- 
actions: in this case the insensitivity of the Higgs mass to scales up to A, become solutions 
to the hierarchy problem, since by the definition of renormalizability, A can be taken to 
infinity. We explore this in sec. |3| and sec. |j. We are led to consider models of the weak 
interactions based on the Hosotani mechanism [34-36]. 

We show however in sec. ^ where the Hosotani mechanism, and the fermionic sector 
are considered in detail, that the simplest such model, based on a bulk SU (3) gauge theory, 
is ruled out. In particular we show that in order to get the right Yukawa interactions, we 
are forced to add sufficiently many bulk matter fields that the SU(3) Yang-Mills is no 
longer asymptotically free in four dimensions and thus has no continuum limit within the e 
expansion of 4 + e dimensions. As we note in our conclusions (sec. ||) this only means that 
one should try to construct more involved models, for example with a larger gauge group, 
and in addition introduce the bulk fields more economically. 

We start in sec. ^ with a short explanation of the e expansion, comparing to its use in its 
more traditional context. We then review the evidence such as it is for higher dimensional 
fixed points. After this we use the known four-loop four-dimensional beta function [37, 38] 
to investigate the e expansion for these fixed points. (It seems that this is the first time 
the early two-loop investigation [23] has been extended.) 

2. Renormalizability in extra dimensions 

Scalar field theory in three dimensions has a Wilson-Fisher fixed point [39]. In the case of 
an iV-component scalar field theory with O(N) invariance, tuning the couplings so that the 
theory lies close to this fixed point results in universal behaviour in a universality class de- 
termined only by N (see e.g. ref. [40]). The 'distance' from the fixed point sets a mass-scale 
m (the inverse correlation length) which is in fact the only parameter left in the universal 
physics. 2 Exactly at the fixed point, one obtains an interacting conformal field theory. 
The properties of the conformal field theory determine everything that is universal about 
this fixed point, for example the scaling dimensions of its operators determine the various 
critical exponents v,u,ri etc. , pure numbers that are amongst the simplest quantities to 
calculate and/or measure, while the effective potential for cf) yields the universal equation 
of state. Viewed in this way, this conformal field theory and its leading relevant deforma- 
tion constitute a concrete example of non-trivial quantum field theory, i.e. a non-trivial 
continuum limit, about a fixed point other than the Gaussian fixed point. 

2 There is one relevant coupling, which can be parametrized by m, and there are no marginal couplings. 
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As we will see, in a closely analogous way, non-trivial continuum limits may well exist 
for Yang-Mills theory in higher than four dimensions. In order to gain a complete apprecia- 
tion of these phenomenona one needs to use the language of the Wilsonian renormalization 
group [39]. For clarity we remind the reader of the basic elements. One works within an in- 
finite dimensional space of bare actions that include all possible local interactions allowed 
by the symmetries of the theory. In this space, there is the so-called critical manifold, 
which consists of all bare actions yielding a given conformally invariant continuum limit 
(for example the O(N) conformal field theories above or the non-perturbative gauge theory 
cases we are about to discuss). Any point on this manifold - i.e. any such bare action - 
flows under a given Wilsonian renormalization group towards its fixed point; local to the 
fixed point, the critical manifold is spanned by an infinite 3 set of irrelevant operators. The 
other directions emanating out of the critical manifold at the fixed point, are spanned by 
relevant and marginally relevant perturbations (couplings with positive and vanishing scal- 
ing dimensions respectively, but where the latter grow as we move to lower energies, i.e. 
are asymptotically free). Next, we choose a bare action with sufficiently many operators of 
the right type to intersect the critical surface, for some choices of couplings. Note that in 
general we have no way of knowing what type of bare action will do: for a non-perturbative 
fixed point we cannot rely on naive dimensional analysis to classify which operators are 
(marginally) relevant. Instead we simply need to search the full space. 

Now in the bare action, we shift a coupling a little bit away from the critical manifold. 
The trajectory of the Wilsonian renormalization group will to begin with, move towards the 
fixed point, but then shoot away along one of the relevant directions towards the so-called 
high temperature fixed point which represents a theory with only infinite mass scales. 

To obtain the continuum limit, and thus a finite mass scale, one must then tune the 
bare action towards the critical manifold and at the same time, reexpress physical quantities 
in renormalized terms appropriate for the diverging correlation length. 

To confirm the existence of such a continuum limit, one really has no choice but to 
follow through the above procedure for example within the framework of a lattice gauge 
theory computation in higher than four dimensions. 

Unfortunately, to date all such investigations have been limited, with ambiguous con- 
clusions. Thus in ref. [24], an exact renormalization group treatment was used, which taken 
at face value implies that the non-trivial fixed point exists for SU{N) Yang-Mills in D = 5 
dimensions for N < 5. The author computed the critical dimension D cr above which the 
fixed point disappears, finding 5 < D cr < 6 in all these cases. Actually, in ref. [24] N = A 
is not reported but the results for N = 2, 3, 5 do suggest that the critical dimension D cr 
behaves smoothly with N as expected. The problem with such a treatment is primarily 
that there is no control of uncertainties due to truncation of the flow equations. There are 
a small number of lattice studies [25-28]. These are subject of course to limitations from 
systematic and statistical errors, which rapidly become more severe with increasing dimen- 
sion (and increasing N). The general conclusion seems to be that a non-trivial continuum 
limit does not exist in D > 4 dimensions for the simple Wilson plaquette action (at least 

3 this is strictly an assumption on our part for the gauge theory cases of interest 
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for N = 2 and D = 5,6). However as we have already emphasised, there is no reason to 
expect that this simplest action is the correct one in this case. Ref. [28] displayed evidence 
for a non-trivial continuum limit even with the simple Wilson plaquette action, providing 
the extra dimensions are small enough, but it is reasonable to assume that this arises from 
lattice artifacts [24,31]. Ref. [26] reported on the addition of an adjoint representation 
Wilson plaquette but the results were inconclusive. Only in ref. [27] were negative conclu- 
sions drawn from more general actions than the simple Wilson plaquette action, but this 
was in the special case of D = 6 dimensions with large N (=27 and 64) after reduction to 
a twisted Eguchi-Kawai model. 

If these non-trivial fixed points exist, one of the most powerful ways of deriving the 
physics resulting from these fixed point behaviours is by computing perturbative quantum 
field theory in four dimensions and using the so-called e expansion [40]. 

The idea is as follows. In four dimensions we know that the perturbative Yang-Mills 
theory is based around the Gaussian fixed point, and parametrized by a small marginally 
relevant coupling g (just as perturbative scalar field theory is parametrized by a small 
marginally irrelevant coupling A and a relevant mass term). All the independent higher 
point interactions (field strength cubed etc. , in the scalar case 6 etc. ) are irrelevant; their 
effect just amounts to a finite renormalization of the low energy couplings. Because the 
theory is free at the associated fixed point, the scaling dimensions of the irrelevant operators 
are equal to their engineering dimensions, and thus are 6 or greater. Now, providing that 
the scaling dimensions of these operators are continuous functions of the dimension D, it 
must be that they are still irrelevant in D = 4 + e dimensions, for small enough e. Thus 
for small enough e, the usual bare action and renormalized effective Lagrangian must still 
yield the right description. 

In the scalar case, this results in the classic argument for understanding the existence 
of the Wilson-Fisher fixed point: the beta function takes the form 

t id li \ = e\ + 0(\), (2.1) 

where the first term is classical and arises simply because we use the Wilsonian coupling, 
namely the dimensionless combination A = A//, while the second -quantum- term can 
coincide with the beta function in four dimensions (as we will see explicitly later). Since 
below four dimensions the classical term is negative while the quantum term starts at 
higher powers than A and is positive (reflecting triviality of four dimensional scalar field 
theory), we get a zero at some point A = A*, at least for sufficiently small e, corresponding 
to an infrared attractive fixed point for A. 

In the gauge theory case, we have qualitatively the same situation but with the signs 
reversed. The quantum term can be chosen to be the four-dimensional beta function and 
is now negative (corresponding to asymptotic freedom) and this can be balanced by going 
above four dimensions, where the classical term is positive (corresponding to a negative 
dimension coupling g). Again, we get a fixed point g = g*, for sufficiently small e, but this 
time it corresponds to the desired ultraviolet attractive fixed point. Indeed we can now 

4 See also ref. [24]. This is independent of a possible multicriticality of the fixed point however. 
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compute this fixed point and its universal consequences by the e expansion, i.e. simply by 
solving for g* order by order in e. 

In fact for the gauge theory, e expansion of non-linear sigma models is even more 
closely analogous [40-42]. O(N) invariant non- linear sigma models in two dimensions are 
also renormalizable and asymptotically free. By the same logic, we expect that in D = 2 + e 
dimensions, there is still a continuum limit for these sigma models but now based around a 
non-trivial ultraviolet attractive fixed point. This fixed point indeed exists and is nothing 
but the infrared description of the O(N) Wilson-Fisher fixed points already discussed (as 
may be confirmed for example by noting that the large N expansion of this fixed point 
coincides with the large N expansion of the one obtained in D dimensional scalar field 
theory) [40]. 

(There is some controversy over the reasons for the poor accuracy in practice of the e 
expansion for these non- linear sigma models [41,42]; some crucial non-perturbative effect 
is perhaps missing. At any rate the qualitative conclusion is correct, namely that ultra- 
violet fixed points exist above the upper-critical dimension where these theories cease to be 
perturbatively renormalizable. It remains to be understood whether the issues of accuracy 
have any bearing on the case of Yang-Mills we discuss here.) 

With respect to the use of the e expansion in all these examples, a word of caution is in 
order about taking the description too literally. Although the e expansion can furnish im- 
pressively accurate answers for universal quantities [40] (the renormalized quantities, which 
are all we are interested in), of course the four dimensional field theory has very different 
properties in general from the D ^ 4 theory (which latter needs to be understood from a 
Wilsonian renormalization group perspective, as already emphasised). Thus, for example 
for the O(N) scalar field theory in three dimensions, perturbation theory about the Gaus- 
sian fixed point generically includes the now marginal </> 6 interaction. More importantly, 
perturbative three dimensional cj) 4 field theory is superrenormalizable. Its ultraviolet di- 
vergences do not give rise to the four dimensional f3 function in ( |2.1| ). Indeed A suffers 
no ultraviolet divergences at all; only m 2 receives divergent corrections at one and two 
loops. Instead, its massless limit is plagued with infrared divergences, a signal that we 
are working about the wrong fixed point. In 4 + e dimensions, with —1 < e < 0, even for 
arbitrarily small e, cp 4 scalar field theory is still superrenormalizable and its massless limit 
suffers infrared divergences, first arising at ~ — 2/e loops (as follows simply from power 
counting) [40]. In gauge theory in five dimensions we cannot a priori exclude any of the 
full infinity of higher dimensional gauge invariant operators from a valid bare action. The 
theory in five dimensions is perturbatively non-renormalizable; it is plagued with infinitely 
many new types of ultraviolet divergence. In particular, its ultraviolet divergences do not 
give rise to the four dimensional (3 function in (|2.4|). Even in 4 + e dimensions, with arbi- 
trarily small e > 0, Yang-Mills theory has non-renormalizable divergences that appear first 
at ~ 2/e loops (again simply by power counting). 

However, the e expansion allows us to exchange the limits (e —* and g — » 0) and 
follow perturbatively the relevant fixed point, and thus in this way access the required 
continuum limit (i.e. the renormalized quantities) directly. 

We now turn to the details. We write the gauge field in D dimensions as Am = ■^-m^i 
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where the T a are generators of some simple gauge group G, with standard normalisation 

tvT a T b = 5 ab /2. (We use a metric of signature H • • •.) Writing the covariant derivative 

as Vm = 9m — iAu-i the held strength is Fmn = i^M, Vjv]- In D = 4 dimensions the 
Lagrangian density takes the usual form 

£ = --L tT F MN F MN . (2.2) 
2g l 

As we have already stated, we expect that the situation is continuous in e and thus in 
D = 4 + e dimensions, fl2.2|) is still the right description for small enough e. At this point 
it is helpful to recall a textbook argument. In dimensional regularisation and minimal 
subtraction, we replace the coupling in ( [2.2| ) by go, and use the dimensionless renormalized 
g = ii e l 2 g to derive 



-e/2 



+ £;M> , (2. 3) 



5o = /i ' ^5 

where the Z n (g) = 0(g 2n+1 ) are the usual power series in g with numerical coefficients 
chosen to cancel all the poles that arise order by order in perturbation theory as e — > 0. 
Now we differentiate ( |2.3| ) with respect to ^i. Since the theory is finite when expressed in 
terms of g, and using the fact that go does not depend on /x, we can invert the result to 
obtain 

»d„g = ^g + /3(g), (2-4) 

where /3(g) = Z±(g)/2 — gZ[(g)/2, and all the Z% with i > 1 are determined in terms of Z\ 
by the requirement that all the 1/e poles cancel in ( p.4[ ). 

As promised, we see that in a little more than four dimensions (e > 0) we can balance 
the first positive term against the negative power series (3(g) to get a fixed point g*, whose 
properties we can now compute order by order in e. 

For what ensues, it is especially important to note that even though (2.4) expresses 
the flow of g in more than four dimensions, (3(g) is precisely the usual four dimensional 
beta function. Thus to deduce the properties of the model in greater than four dimensions, 
we only need renormalize the model in the normal way in four dimensions and then use 
Q2.4D and similar equations as a formal device to analytically continue to e > 0. 

As we will now show, the e expansion for these ultraviolet fixed points is in fact 
surprisingly well behaved. We can illustrate this using the known four-loop (3 function 
[37, 38] for SU(N) Yang-Mills. 

At this point we note a convenient point of detail: the four-dimensional (3 function is 
the same in the MS and MS schemes [38]. This can be understood from the definition of 
MS namely that, instead of subtracting via counterterms, powers of 2g 2 /(4-/r) 2 times poles 
in e, one substracts powers of 

_^!_ {1 + ln(47r ) e - 7£e} 

(jE being Euler's constant), the corrections arising from 

2 



(4^)0/27^/2)' 



-7- 



the D dimensional solid angle divided by (2tt) d that accompanies each loop integral. There- 
fore for the Z n in (|2,3j), MS amounts to an e dependent rescaling of g. When working at 
finite e, the results can be expected to be better behaved if we instead absorb the full 
dependence coming from [40]. This amounts to replacing 2g 2 /(47r) 2 in the MS scheme, 
with g 2 0,£). Clearly this modified MS also yields the same four-dimensional (3(g)- Finally 
for convenience in the ensuing analysis, we also absorb a factor of N and write 

a:= Nn D g 2 = NQ D n e g 2 (2.5) 

In terms of a, ( |2.4D becomes 

fid^a = ea + (3(a), (2-6) 

where 

/3(d) = -(3 a 2 - (3 X 6? - (3 2 a 4 - (3 3 a 5 + 0(d 6 ); (2.7) 

these (3 n are the coefficients quoted in ref. [38] but divided by 2 n N n+1 . 

We compute a*, although this is not universal, but we also compute the index v via 

[23] 

;==-!«•> =-!<*•>• <«> 

which is the universal power relating the analogous QCD-like scale Aym to the distance 
from the fixed point at any high scale \x where the coupling is a(/j,) is sufficiently close to 
d*: 

Aym oc fi\a — a*\ u . (2.9) 

Solving (|2.6|) for 5*, 



(3o (3$ 



[ 2 PI (3t) e + { 5 (31 (31 5 (3l) e + ° (e) ' (2 - 10) 



d, = ^-^+(2^-^U a +(5 



and substituting the values for (3 n , we find a surprise. Whereas the coefficients (3 n show the 
expected dramatic increase in n, consistent with ( |2.6| ) being asymptotic (even at N = oo): 

/3b = 3.667, (3t = 5.667, ^2 = 13.23, (3 3 = 39.43 + 51. 22/iV 2 , (2.11) 

the series for d* 

d* = .2727e - .1150e 2 + .02372e 3 - (.007395 + .07729/iV 2 )e 4 , (2.12) 
is considerably better behaved, and at N = oo it is not at first clear that this series is 



asymptotic. Substituting ( 2.10 ) in fl2.8j), we obtain 



- = e + .4215e 2 + .1813e 3 + (.1242 + .8502/iV 2 )e 4 . (2.13) 

Again, at N = oo, this series looks very well behaved. 

We start by investigating this limit in more detail. Writing the coefficients of e n in 
fl2T2D and (glaD and v n respectively, successive ratios |a n /a n +i| give 2.37,4.84,3.21 
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for n = 1, 2, 3. And similarly successive ratios \v n / l 'n+i\ g iVe 2.37, 2.32, 1.45. If we did not 
know better, we might even be tempted to conclude that these ratios have a limit r > 0, 
corresponding to a finite radius of convergence r. However, we see from (2.12), that d* 
has negative values for all e < 0, as expected, since the four dimensional beta function 
analytically continued to a < behaves like that of a trivial theory and thus allows us 
to balance the two terms in (2.6) and find a fixed point in D < 4 dimensions. Since the 
theory is non-perturbatively sick for a < [43], the above series in e must have zero radius 
of convergence. 

If the series are already displaying their limiting asymptotic behaviour (viz. a n ~ 
cs n T(n + Q/R n , for s = —1 and some constants c, ( and R, and similarly for u n , with 
s = l) 5 then we can estimate them by taking the sum to the point where the terms stop 
getting smaller (including the smallest term) and using the next term to estimate the error. 
In this case we ought to see \a n /a n+ \\ ~ R/(n + 0, and thus lrm n ,_» 00 n\a n / a n+ \\ = R. 
Forming this combination for n = 1,2,3, we get for the a n s: 2.37,9.69,9.62, and for the 
v n s: 2.37,4.65,4.38. 

We can already draw a number of important conclusions. Firstly, there is clear evi- 
dence that the series do approximate their limiting behaviour already after the first term. 
Secondly, the fact that the last two combinations are so close to each other, for both series, 
indicates that the corresponding Q must be small in these cases; of course the series are 
not long enough to reliably estimate it however. Thirdly, we expect the ~ i? th coefficient 
to be the smallest one, thus for \jv in ( |2.13| ), we expect the next (5-loop) term to have a 
larger coefficient than the four- loop term, while for d*, we expect to have to go to 10 loops 
in (2.12) to see the coefficients start growing with n. 

Fourthly, we can use the series to estimate v and d* for positive integer e, as follows. 
For v we simply estimate it as described above, from the asymptotic series for \jv\ for 
e = 1 we have a special case in that we are missing the 5-loop term that would provide an 
estimate of the error. This is why there is no error for the appropriate entry in table ^. 

For d* at e = 1 we have an alternating series in which we expect the terms to keep 
getting smaller until the (R — l) th term with the slightly larger R th term providing the 
error, where R ~ 10. Clearly in this case d* should lie between the sum of the first three 
and the sum of the first four terms. 6 This gives us the D = 5, N = oo estimate in table 
|l[ In general the point where the terms stop getting smaller can be expected to be at 
n = R/e. Comparing with (2.12), we see that for integer e > 1 we should use the above 
method for summing asymptotic series, except that for e = 2 and 3, we are missing the 
term that would provide the error. This is why we do not display an error when D = 6 
and N = oo. 

Finally, at finite N, the series are clearly asymptotic, the 1/N 2 terms in ( |2.12 , 2~T3| ) 
(the result of the first divergent non-planar contribution - which appears at four loops [38]) 
having a much larger coefficient (and of the same sign). Without further information (see 
below), we can estimate the sums only in the cases where this last term is not the smallest 
term. We find that the last term is smallest in ( |2.13| ) only when D = 5 and N > 4, and in 

5 equivalently a n ~ c(—) n n c n\/R n 

e If the coefficients kept getting smaller this would be a theorem. 
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D\N 


2 


3 


4 


5 


oo 


5 


.18(3) 








.178(4) 


6 


.28(43) 


.28(26) 


.28(20) 




.157 



Table 1: Estimates for a* for varying N and D. 



J2.12 ) when D = 5 and N > 3, and when D = 6 and TV > 5. This is why we the leave the 
corresponding entries blank in the tables. 

In addition, we find that for 
l/i/, the terms in the series only 
increase once D > 7 (for any N). 
For a» when D = 7, the error is 
larger than the estimate for 2 < 
N < 12, after which the last term 

is the smallest so that we cannot provide an estimate. At N = oo, we saw above that we 
cannot estimate the error, and the series sums to d* = —0.175, which cannot make sense 
physically. For D = 8, the error is much larger than the estimate for d* for all N. For 
D > 9, the terms only increase in magnitude with n. We take all this as evidence from the 
e expansion that the fixed points do not exist in D > 7. 
In D = 6 the evidence is some- 



D\N 


2 


3 


4 


oo 


5 


.62(13) 


.62(9) 




.579 


6 


.19(20) 


.19(13) 


.19(11) 


.19(8) 



Table 2: Estimates for v for varying N and D. 
oo the error is 40%. For d* the errors remain large where they can be estimated 



what marginal. At N = 2 the error 
is larger than the estimate, both for v 
and d*. Perhaps this indicates that the 
fixed point does not exist. For v the er- 
rors are large but gradually decreasing 
for N > 3 and N increasing, but even 
at N = 
at all. 

On the other hand, for D = 5 we appear to have clear evidence for a fixed point: the 
errors are small where they can be estimated. There are missing entries only because the 
series are too short in these cases. Given the other values it is reasonable to assume that 
d* ~ .18 and v ~ .6 for all N > 2. 

These conclusions are in broad agreement with the earlier alternative approaches al- 
ready discussed [24-28], given the limitations with all these studies. 

We have only performed the simplest estimates. It should be possible to do much 
better. In the lower dimensional scalar field theories, powerful techniques have been devel- 
oped to cope with their divergent e expansions. One first transforms the series to the Borel 
plane. Relying on no other knowledge the series can be resummed by Pade approximants. 
However, if one assumes the e expansion inherits the large order behaviour of A</> 4 theory in 
three dimensions one knows about the closest singularity in the Borel plane (from studying 
Lipatov instantons) [40]. By conformal mapping the Borel plane, this information can be 
taken into account. The results are impressive: for better known quantities e.g. A* and 
v (and using also information from the exactly soluble two-dimensional Ising model) the 
theoretical error has been reduced to a few per mille [40] . 

In the gauge theory case, one can also transform the series to the Borel plane and 
use Pade approximants. However, one also has similar knowledge about the large order 
behaviour of perturbation theory coming from the closest singularities in the Borel plane, in 
this case due to infrared renormalons and instantons, which can also be taken into account. 
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3. Renormalizability with matter and branes 

We pursue theories which have non-perturbative continuum limits in more than four space- 
time dimensions, in the sense of having ultraviolet fixed points within the e expansion. We 
have seen that this requires that when e = 0, the model must be renormalizable and in 
addition all couplings in the bulk are relevant (effectively masses and couplings of positive 
mass dimension) or marginally relevant (equivalently asymptotically free). 7 

Thus we can readily incorporate non-Abelian gauge fields in the bulk of the higher 
dimensions. We can then add fermions to the bulk of the higher dimensions, providing 
there are not too many or that their representations are too large, such as would destroy 
the coupling's asymptotic freedom. It is harder to include scalars in the bulk (because the 
scalar self-coupling must also be asymptotically free), but possible for a careful choice of 
couplings. On the other hand, we cannot at the same time as keeping asymptotic freedom, 
use the scalars to give a mass to all the gauge fields via the Higgs mechanism [44] . Inclusion 
of both fermions and scalars allows asymptotic freedom only for a fine tuning of the Yukawa 
couplings which is natural only in supersymmetric theories [44]. Finally, Abelian gauge 
fields can be included only on four dimensional branes if at all. 

As is already clear, these models are highly restricted. We will however meet yet other 
constraints renormalizability puts on the couplings. We will concentrate on the case of five 
dimensional theories, i.e. where one should set e = 1 eventually. This is because, as we 
have seen, the evidence for non-trivial fixed points is strongest in five dimensions. However, 
most of our comments obviously extend to other dimensions. Our conventions are to write 
five dimensional indices in Latin capitals with components M = 0, 1, 2, 3, 5; coordinates we 
will write as x M = while we write momenta pm = ip^,Ph), the bar being an extra 

flag for four-dimensional parts. We write the Lagrangian density of the non-Abelian gauge 
field as in eqn. fl2.2|) . 

We do not consider gravity since we do not know how to describe it as a renormalizable 
theory. We thus take the conservative view that there is a flat classical gravitational 
background, to be set and not questioned. Phenomenologically therefore, these models 
have to have the fifth dimension compactified, which we can take to be a circle radius R: 
y is identified with y + 2irR. Although the models we consider, are valid as quantum field 
theories up to infinite energies, in practice we have to set the cutoff at the natural five- 
dimensional Planck mass since quantum gravitational effects cannot be ignored above this. 
(This means that A ~ a 1 / 3 x 10 13 GeV, where we have taken the radius of compactification 
to be 1/R = a TeV [45].) 

At first sight it is attractive when computing such a model in the e expansion, to 
regard the extra e dimensions as the ones that are compactified. This is not consistent 
however. In the e expansion the extra dimensions only enter by the change in the scaling 
of the classical coupling constant. The calculation is thus not sensitive to topological or 
geometric features of the e extra dimensions. This is reasonable: the approach to the 
fixed point g*, is a property of the theory in the far ultraviolet, i.e. at very small length 

7 Note that these are the same conditions that are believed to be required to ensure the non-perturbative 
existence of the quantum field theory in four dimensions. 



- 11 - 



scales, where space-time should look flat. Where a more precise matching is required 
between renormalized couplings above and below a compactification scale, the result from 
e expansion with one of the original dimensions compactified, should clearly be used to 
compute threshold corrections. 

If parts of the five dimensional gauge fields Am are to represent some Standard Model 
gauge fields, then we need to go beyond trivial compactifications. (We do not observe the 
adjoint massless scalars A5 that would arise from such compactified Yang-Mills fields, and 
we cannot use extra-dimensional Higgs fields to give masses to all parts of ^,5.) There are 
two overlapping possibilities that have been suggested in the literature, the extra dimen- 
sional space can be an orbifold or more generally can have a boundary on which non-trivial 
boundary conditions are imposed [14]. 

Either way, once we take the model seriously as describing a renormalizable continuum 
limit, we have to introduce boundaries, a.k.a. branes, and boundary couplings which we 
cannot ignore, as we now explain. 

In the general case [14], boundary conditions on the gauge fields are imposed that 
explicitly break the gauge group G down to some smaller {e.g. Standard Model) group 
H. Of course an explicitly broken gauge theory is not renormalizable. This is the only 
problem, and is not in fact a problem per se, for the little hierarchy solution envisaged by 
these authors, as can be understood by viewing the boundary conditions as imposed by 
boundary Higgs fields in the non-linear sigma model limit [46] . The effective cutoff for the 
sigma model can be identified with the effective cutoff the theory already necessarily has. 
However these problems do rule out general breaking by boundary conditions as an option 
here, because the non-linear sigma model is not renormalizable. 

The orbifold alternative consists in dividing out the extra dimension by a discrete sym- 
metry which is not freely acting. For illustration we consider only the simplest case of Z2 
parity y h-> —y. 8 We are interested in the case where the parity does not commute with the 
gauge group (otherwise all of A§ will be odd, thus gain Kaluza-Klein (KK) masses of order 
the compactification scale, and just return us at low energies to a four dimensional world 
of Yang-Mills with gauge group G). The gauge fields satisfy Am(x,v) = PAm(x, —y)P- 
(It is to be understood that P also maps any 5th Lorentz component to minus itself.) Pro- 
viding P generates an automorphism of the Lie algebra of G, such an action is consistent, 
but the result is that the gauge symmetry is typically restricted to a smaller semi-simple 
group H at the orbifold points y = 0, irR. Indeed, we can always choose a basis for the Lie 
algebra so that the action of P is diagonal. The even generators, PT a P = T a , generate the 
subgroup H, whilst all the other generators T a ' are odd. Although from the low energy 
four-dimensional point of view, the gauge group G appears to have been broken to H, this 
is not really true. Gauge transformations 5Am = Vm^ are only restricted in the sense 
that Q must also satisfy Q(x,y) = PQ(x,—y)P, and thus the components Q. a ' vanish at 
the boundaries y = 0, irR. 

As noted in ref. [47] , the orbifolding of the 5th dimension results in divergent quantum 
corrections which are localised at the orbifold points. We are thus forced to introduce brane 

8 Most of our comments apply equally well to more involved cases, the most general case in this situation 
being Z 2 xZj. 
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Figure 1: One-loop contribution to the (3 func- 
tion in background field gauge. The internal line 
stands for all components, including ghosts. 



Lagrangians at these points even if we did not start out with them, with couplings of the 
same form as the divergences (so that renormalizing the couplings absorbs the divergences). 
By the symmetry at the branes, we need four times as many couplings 1/gf as there are 
semi-simple factors in H, to multiply separately the generated terms (J~^ u ) 2 and (J 7 ^) 2 
on each boundary (see sec. Q for an example). Using the Lagrangian in the form (|2.2|), 
and for example working in background field gauge [48] cf. fig. [[], it is easy to see that the 
couplings then run according to ^8^(1 /gf) = 6j/(87r 2 ) where these bi are combinations of 
group theory factors, yielding numbers say, in the range 0(1) to 0(10), but whose value 
depends on the details (in particular the orbifolding and the fermions we are about to 
introduce). 

Note that in computing these effects 
we work in 5 dimensions: the localised di- 
vergences only arise from integrating over 
the four-dimensional loop momentum ft**, 
whilst the 5th component only participates 
in the Umklapp process ~ $2k 5 ,p 5 ±p' s which 
in position space leads to the brane localised 
delta functions 6(y) + 5(y - ttR) [47]. 9 We 
would get the wrong answer if we worked in 
4 + e dimensions, thus expanding the k mo- 
mentum integral around three dimensions. 
Indeed in this case the diagrams would all 
be ultraviolet finite, and all infrared finite, 
except for the k§ = term which is linearly IR divergent for vanishing p. Whether we work 
in 5 dimensions or expand around 4 dimensions has to be determined by the effect being 
calculated: we need to choose to expand about 4, if and only if it is the critical dimension 
for the effect in question, i.e. the dimension at which logarithmic divergences appear. At 
one loop, this clearly decides the issue. Perhaps at two loops and higher, there can be an 
effect that can only be properly understood by considering its divergences both in 5 and 
around 4 dimensions? We leave this question for the future. 

Note that it is not possible to argue that the 1/gf are much less than bi/(8n 2 ), because 
the result would then be strongly [i dependent. Indeed in this case, changing the value of 
H by a factor of two, would return the couplings to l/gf(p) ~ 6j/(87r 2 ). We thus conclude 
that the brane couplings are O(10 _1 ) or greater. 

Even if we were to limit the appearance of these brane kinetic terms by arranging the 
bi to be very small or cancel (at one loop or to any number of loops), it is not possible to 
set 1/gf to zero if there are other non- vanishing brane interactions involving A, without 
generating further divergences, as we explain below. 

Since these couplings 1/gf do not respect the full gauge symmetry G, Grand Unification 
in a single gauge group is not really possible in these scenarios (even if some of its features 
such as charge quantization can be preserved [18]). For this reason, and the fact that lack 
of renormalizability is less of a problem in any case at Grand Unified scales (cf. sec. |l|), we 
do not pursue the possibility of reno rmalizable extra dimensional Grand Unified Theories 
9 The loss/gain of KK momentum together with possible reflection of the momenta is analogous to 
Umklapp interactions with lattices in condensed matter. Indeed in reality the branes absorb the change in 
momentum. 
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further here. 

To see that it is not possible to switch 
off brane kinetic terms for A in the presence 
of other non-vanishing brane interactions for 
A, consider adding an H invariant scalar field 
4>{x) to the brane(s). This will generate from 
the kinetic term in particular a brane interac- 
tion of the standard "seagull" type 




Figure 2: Contribution to branc-localised 
5{y)AA4>4> (3-1) interaction. The dashed lines are external 0s 

The dotted lines are insertions of external y- 



4 



2 - J (3.3) 



to be added to the overall Lagrangian density momentum from (|3.2|). 

J2.2| ). The effect we are about to describe happens for any field propagating in the bulk and 
any interaction, so for the moment we will ignore the Lorentz and colour indices, treating 
A as though it were a scalar field. Since this is a small-distance effect, it occurs just as 
well in infinite flat bulk dimensions with a single brane at y = 0. 

Now the one-loop diagram indicated in fig. § generates a local divergent correction to 
the brane ^(x) interaction. The delta function in ( |3.l[) is easily taken into account if we 
regard it as a background field 

Mx,v) = 8(y) (3-2) 

in a five-point interaction, as indicated in fig. [| We see that the local interaction generated 
is actually 

H 

which is not well defined with the identification (|3 

Part of the problem of course lies in our assumption that the brane is infinitely thin 
and infinitely heavy. In a fully realistic situation, the brane would have a finite mass which 
we identify with the scale of quantum gravity A, and a form factor 

<kfav) = Mv), (3.4) 

where becomes a delta function only in the limit A — > oo. If A were not a gauge field, 
we could indeed replace 6(y) by ([O]). Since for functions smooth at distances 1/A, 

4>l = ckfo, (3.5) 

(where c is some constant of order 1) this extra divergence could be absorbed into the 
couplings. However, in our gauge theory, the replacement of 5{y) by a form factor would 
be fatal: in order for ( |3.1| ) to remain gauge invariant, (ft has to be extended so that it 
depends on the bulk dimensions. This would mean we would not have the option to 
have brane-localised fields, severely limiting not only the model-building possibilities but 
also the possibility to keep the bulk theory asymptotically free in four dimensions and thus 
supportable by an ultraviolet fixed point in 4+e dimensions. An alternative way to proceed 
might be to use ([^) and promote relations such as (|3.5| ) to identities. At this stage it is 
not clear if the model can then be made fully renormalizable. 
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Fortunately we do not have to pursue these directions further, since the divergence 
does not appear if we keep a brane kinetic term for A. To see this, again for simplicity we 
treat A as a single component real scalar field. After a rescaling the kinetic terms yield 
the inverse propagator 

A" 1 = - d M d M - - 5(y) 8^, (3.6) 
m 

where m = 2gf /g 2 is the ratio of couplings and has dimensions of mass. The corresponding 
propagator has already been derived in the literature (see e.g. [49]). However, to discuss 
divergences it is helpful to write it completely in momentum space. This is easily found 
directly as follows. In Euclidean momentum space, 

P 2 

A~ 1 (p,p 5 ;<L<75) =p 2 S(p-q) + 2 — 5(p-q), (3.7) 

m 

where 8(k) (5(k)) means the D = 5(4) dimensional delta function multiplied by (2tt) d . 
Substituting the ansatz, 

A(p,J> 5 ;<?,<?5) = a(p)S(p - q) +b(p,p 5 ,q 5 )5(p - q), (3.8) 

where a(p) and b(p,p5,qs) = b{p,q^,p^) are functions to be determined one finds 

A(p,p 5 ;q,q5) = 2 21-2 \ 2\t TFh- ( 3 ' 9 ) 

(The modulus term \p\ = s/p 2 arises from integration over 55.) 

Taking account of momentum conservation in the usual way, one finds the contribution 
fig. H is proportional to 



d 4 xdV^(x,0)^(x',0) J-J^ e W- x > j — ^A(q)A(q + p), (3.10) 

where the effective four-dimensional propagator 

A(p) = 1 (3.11) 
p z + m\p\ 

comes from integrating over the injection of r§ (r' 5 ) with coefficient 1, that arises at each 
vertex from ( |3.2| ). 

Note that the effective propagator results in no worse ultraviolet divergences than a 
four dimensional theory completely confined to the brane. It is clear that all one-loop 
diagrams of this sort will be regulated in the same way and thus no terms of the form 
fl3.3j ) are generated. Intuitively one can understand what is happening as follows. For very 
small distances on the brane, the S(y) term in ( |3.6| ) dominates. In this regime everything 
is confined to the brane and we have just a four dimensional theory. The five dimensional 
coupling constant g, through m, provides a cross-over scale in the effective propagator 
( |3.11 ), so that only for distances larger than 1/m do we feel the effects of the full 5 dimen- 



sions. In fact what we have is the gauge theory analogue of the "DGP" model [49], where 
these effects are discussed in a gravitational context, a neat derivation of precisely ( 3.1 1| ) 



- 15 - 



appearing in ref . [50] . This intuition gives us confidence that this ultraviolet regularisation 
effect works to all loop orders. 

Finally, note that the above considerations imply that all the functions of the brane- 
localised fields (here that of a brane scalar field) must be singular in the limit 1/gf — > 0, 
in which the brane kinetic terms are turned off. 

Note that the gauge field components Ac- are even under P. Their lowest KK modes 
thus behave like scalars as far as the low energy four-dimensional theory is concerned. 
There is nothing in the four-dimensional theory that prevents them from gaining a mass. 
From the five dimensional point of view this is not allowed by gauge invariance. However, 
an effective action can and will appear for the non-local operator 

§ p (x) = $(x,2ttR), (3.12) 

the Wilson line that winds around the compact dimension, where 

§(x,y) =Pexp -i [ V dyA 5 (x,y) (3.13) 
J o 

is the Wilson line from y = to y = y, located at x (and V stands for path ordering). 
Since under gauge transformations, 

5$(x, y) = i Q(x, 0) $(x, y) - i $(x, y) fl(x, y), (3-14) 

3>p(x) transforms homogeneously as an adjoint scalar at (x,0). 

Since $p contains a power series in A§ , the gauge invariant effective potential trV e ff(<frp) 
not only generically results in masses for the lowest KK mode of the components A§ , but 
can also result in spontaneous symmetry breaking: <A§ > 7^ 0. This is the Hosotani mech- 
anism [34]. In the spontaneously broken phase, although we can by a gauge transformation 
set A5 = 0, this is at the expense of a non-trivial Scherk-Schwarz twist in periodicity con- 
ditions for A [51] (as we review in sec. |I|) which also breaks the gauge group. Meanwhile, 
V e ff and the low energy physics as described through it, are completely unaffected by this 
gauge transformation. Note that the mass and more generally the effective potential is 
protected from divergent corrections: the effect is non-local in the full five-dimensional 
space and thus quantum corrections are cutoff by the compactification scale l/(2irR). 

Clearly an extremely attractive possibility now arises, namely to regard the Higgs 
scalar field as the components A§ [6-13,34-36], for if we can bring such a model into 
accord with phenomenological constraints, whilst keeping it renormalizable, we will have 
solved the hierarchy problem. This is the possibility we will consider in more detail in 
sec. |j. We will use the framework developed there to make general comments on model 
building constraints that arise when we consider couplings to fermions. 

4. The simplest gauge-Higgs unification model 

The simplest possibility for the unified gauge group is G = SU(3) [52]. 10 With the gen- 
erators T a = A a /2, where the A a are the Gell-Mann matrices, we can decompose Am 

10 This is really SU(3) x SU(3). The second SU(3) is for strong interactions and is readily incorporated 
in the bulk. Since the real problems lie with the weak interactions we will not discuss colour further. 
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as 



Here, (Wm)} 5 with i,j = 1,2, contracted into T a (a = 1,2,3), are the gauge bosons 
associated to the top left SU(2) subgroup. This will be identified with the SU(2)w part 
of the Standard Model gauge group. Similarly will be identified with the U(l)y gauge 
boson, while K % = H\ = \f2{A<$)\ is to be identified with the SU(2)w doublet Higgs. We 
implement the restriction of SU(3) to SU(2)\y x U(l)y, by identifying the parity operator 
P as mapping the third component of the SU(3) fundamental representation to minus 
itself, whilst leaving components i = 1,2 alone. Under this map indeed, W^, and H 
are even and thus at this stage will have massless KK modes, whilst W§, B§ and are 
odd and thus only have KK masses of order the compactification scale or greater. 

The radiatively generated potential for <I>p, will then spontaneously break the model, 
effectively giving H a vacuum expectation value, via the Hosotani mechanism as outlined 
in the previous section. Using the remaining global symmetry we can ensure the vacuum 
expectation value takes the form < H >= (0, v) / v2 for some real v. Thus SU (2)w x U(l)y 
will be broken down to electromagnetism, just as happens in the Standard Model. 

Since 

[T S ,H] = ^H, (4.2) 

we already have a problem however. If we were to put the coupling g back in its normal 
place in the covariant derivative, we would want to identify gT 8 with g'Y/2 where Y is 
the hypercharge and g' the associated coupling. Comparing to the Standard Model Higgs 
for which Y = 1, we get from (|4.2j), that g' = V3g and thus we 'predict' sin 2 6\y = 
g' 2 /(g 2 + g' 2 ) = 3/4, a phenomenological disaster. 

However, we now recall that we necessarily have the following couplings on the branes: 

AC = 5(y)£ 1 + 5(y-irR)£ 2 , (4.3) 

where, writing out the non- vanishing parts of J- 2 5 and similarly the semi-simple components 



Here we have defined H' = d^H^, the brane couplings w a , b a and h a , and the field strengths 
B^v = d^Bv - d u B^ and = i[D^, D v ), where = 8^ — iW^. Note that 

V^H = d^H - iW^H - i^-B^H, (4.5) 

as a consequence of fl4.2| ) and the restriction to the branes. 

Recall that the couplings 1/w 2 , 1/b 2 and 1/h 2 are 0(1/10) or larger. (The symmetric 
point where they have the same value on each brane will typically be broken when fermions 
are included.) In the presence of such couplings (|4.4|), spontaneous symmetry breaking 
deforms the lowest mass modes for the weak vector bosons and so that they are 
no longer simply constants in the y direction (in any gauge) [13]. 
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However, this effect is controlled by the ratio of scales 8/2n = Rv/2 [13], vanishing as 
6^0. This ratio appears in 



<$ P >= cos6» -ism9 . (4.6) 




The physics is invariant under 6 \— > + 2tt. Defining 8 to be in the fundamental domain 
[0, 2tt), the natural theoretical expectation is that 6/2tt ~ 0.5 if spontaneous symmetry 
breaking takes place. However, indirect limits require the compactification scale 1/R = 
aTeV where a > 2 — 5 or greater [53]. Identifying v for the moment with the Standard 
Model Higgs' vacuum expectation value, we require 9/2ir < 0.06. In this case we can ignore 
the deformations in the first approximation and take the lowest mass modes for W^, 
and H to be constant in the y direction. 

The form of the effective Lagrangian for these modes is then 



2-kR 



' "- ; ' -y+ , 1 IV7 ZT|2 



2g 2 » v 4g 2 ' <7 2 1 ^ 1 

¥ w l» E i - -A E w + 1 v ^" 2 E I . (4-7) 



where ( |4.5D holds but now with y-constant modes. Thus clearly if we define 

1111 , ln , 
— = + — + — (4-8) 

3 111 

~g\ ~ 2^ + q + q { ' 

where we have used the Wilsonian dimensionless bulk coupling ([^) evaluated at [i = 1 /R 
(and D = 5, N = 3), du. a = g 2 /4ir 3 R, and redefine 

s, ; „ ,, / i i . O - 1 ' 2 



= B » = Tz B - H = H {^, + lq + M) ■ <410) 

to put couplings back in their usual places, we get back precisely the relevant part of the 
Standard Model Lagrangian 

-\trW 2 ^-^Bl u + \V,H\ 2 (4.11) 

where in particular W^ u is now defined as expected in terms of = 8^ — igiW^ and the 
full covariant derivative 

V lt = d lt -ig 2 W^-ig 1 -B lt (4.12) 

defines a hyper charge 

Y = diag(l,l,-2)/3 (4.13) 
which gives the Higgs Y = 1 as expected. 
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Taking into account the rescaling to H, we see that v is actually defined in terms of 
the Standard Model Higgs' expectation value v by 



R ~ v =^i^ + h l + h 2 (414) 



These equations allow us to make several straightforward but important conclusions. 
Firstly, in the limit that is sufficiently small to neglect deformations of the lowest KK 
wavefunctions, we regain from ( 4.11[) , the custodial symmetry of the Standard Model and 



thus determine the p parameter to be one at tree level, as required. Secondly, using the 
experimentally determined numbers [54], the bound 1/R = a TeV and the fact that l/w 2 , 
1/b 2 and 1/h 2 are > 1/10 we can place further bounds on the values of the parameters. 

Thus from ( |4.8|) , l/2ir 2 a < 2.3, and with natural values for w a , we have a > 0.02. 
This is clearly perturbative and an order of magnitude smaller than the typical fixed point 
value (cf. sec. ||). At energies higher than 1/R, a will run according to (|2.6|), reaching the 
non-perturbative physics associated with the fixed point at energy scales ~ 10/ R. 

On the other hand we clearly have an upper bound on l/w 2 + l/w 2 of 2.3. From 
Q4.9| ) the strict upper bound on l/2-7r 2 d implies a lower bound on 1/b 2 + l/b\ > 21. We 
therefore require very large U{l)y kinetic terms on the boundary. Perhaps these arise 
naturally from non-perturbative physics close to the ultraviolet fixed point a*, bearing in 
mind that Abelian gauge fields would be separately non-renormalizable in the bulk. Of 
course phenomenologically, the origin of the large values is the discrepancy between the 
'bulk' sin 2 9w = 3/4 and the experimental one. 



The Lagrangian given by (2.2) and (4.4) is very similar to the bosonic sector of the 



Lagrangian considered in ref. [13], except that Scrucca et al. introduce an extra bulk U(l) 
gauge field to allow sin 2 9\y to be set to its experimental value (and do not introduce the 
h a terms). We do not have the option of including a bulk U(l) gauge field since it is not 
renormalizable within the e expansion. 

Finally, (4.14) implies that 9/2n > 0.08/a with natural values for l//i 2 . Therefore 



phenomenologically preferred values of 9 are consistent with our approximation. 

To allow for larger 9, we have to increase a. Thus a ~ 0.1 (implying l/w 2 ~ 1) allows 
for 9/2tt ~ 0.18/a, but to get the theoretically natural values of 9/2ir ~ 0.3 [13] requires 
a ~ 1.2. Such a large value would probably (depending on matter content) imply that the 
fixed point a is being approached from the right (i.e. from the region a > a*). There is no 
problem of principle with this happening but in particular, large higher order corrections in 
the e expansion mean that we could then no longer trust these simple formulae. However, 
in any case we are back in a regime where significant distortions from Standard Model 
relations will be found. We should note that there are not enough parameters in the model 
in this regime to tune away the resulting anomalous triple gauge couplings [13,55], tune 
p, sin 2 9w and mz and the effective brane couplings to matter fields [13] to their correct 
values. Even if we could manage this, it would be an unsatisfactory accident since custodial 
symmetry has been badly broken. 

In summary, for the bosonic sector of the model to be phenomenologically acceptable 
we must have large values of 1/b 2 ~ 10, and a much smaller 9 than we would find without 
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some special mechanism. Indeed we can choose a small a > 0.02 and natural values of 
l/w 2 ~ l/h 2 ~ 1/10, in which case we have to arrange the model to dynamically determine 
0/2vr~ 0.04 -0.09. 

Before turning to the introduction of fermions we discuss briefly the peculiar remnants 
of gauge invariance operating on the brane. Decomposing a gauge transformation similarly 
to flPD, as 



we have that u and (3 are P even, while <p is P odd. Thus u>, [3 and (f)' = d^(j) survive on the 
branes whilst d§uj, d§(3 and <j> vanish on each brane. It follows that the brane Lagrangians 



(4.4) are invariant under the surviving, or "remnant" [11], gauge symmetries: 



SWp = [D^u] 
SB^ = d„p 

/3 

5H = <f>' + iuoH + i^-/3H 

5H^ = iuiH^ + i^-m'v + V M 0' ; (4.16) 

where V^0' is ( |4.5| ) with H replaced by 4>'. 

As noted in ref. [11], the shift symmetry ft protects against brane mass terms appear- 
ing for the Higgs. It is interesting also to note that if it were not for the bulk Lagrangian 



2T4 ), the brane Higgs kinetic term in (4.4) would be trivial since it could be gauged away by 



a finite <f>' = H transformation in ( 4.16| ), H' playing the role of an auxiliary field. However, 



this term becomes non-trivial when considered as part of the full Lagrangian: although we 
can still gauge the brane Higgs kinetic terms away, this is a gauge choice which in general 
conflicts with the need to make other gauge choices (for example the background Feynman 
gauge used typically to compute the radiative potential). Furthermore, H' is of course no 
longer an auxiliary field but part of the bulk degrees of freedom evaluated at the orbifold 
points. 

In order to complete a realistic description we need to introduce quarks, and leptons. 
For example consider initially, the top and bottom quarks. We write their SU(2) W x U(l)y 
2i/6> I2/3) I-1/3 representations as Ql = (ti,6i), £r and 6r. As usual in gauge-Higgs 
unification models, we cannot introduce these as bulk fermions because the interactions 
would be flavour-symmetric and have the wrong hypercharges, following from ( |4.13 ). If 



we introduce them as fields that live only in the brane(s), then we can simply assign them 
the correct hypercharges as according to ( 4.12[ ), since only the restricted symmetry ( 4.16j ) 



is active there. However, ( [4.16 ) as well as including the standard SU(2)w x U(l)y gauge 



transformations (after use of ( 4.10 ) and defining = v3/9) includes the shift symmetry gen- 
erated by (/)' which forbids the Yukawa interactions, for example under the shift symmetry 
we have 5(Q L Hb R ) = Q^'bR. 

(We can cancel this by postulating heavy mirror fields, e.g. a 2jy 6 Qr, adding appro- 
priate mass terms MQlQr (plus c.c), and defining 5Qr ~ (p'bR, and so on, resulting in 
a see-saw mechanism with the lightest states to be identified with the quarks. However 
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we then need kinetic terms for the Qr. The terms generated by the shift transformations 
acting on this can be cancelled by introducing an appropriate interaction with H' and 
appropriate transformations into heavy partners for ba and Unfortunately then more 
fields are needed to in order to cancel new violations of the shift symmetry and so on. It 
appears that it is not possible to find a closure of the symmetry which is non-trivial, linear 
and finite dimensional. However, since we are already introducing an infinite number of 
fields through the KK excitations, it could be worthwhile to pursue the possibility of infinite 
dimensional representations of the shift symmetry. Note that if these could be constructed, 
the shift symmetry would still protect the Higgs from gaining a brane potential, but allow 
standard Yukawa interactions for the real quarks and leptons.) 

This problem has been circumvented in the literature by taking instead the Wilson 
line <£, using <3?p for coupling doublets and singlets on the brane at y = (analogously 
a Wilson line that wraps once round the compact dimension but starting at y = ttR, for 
coupling both representations at y = ttR), or <j?(7ri?) for coupling fermions on one brane to 



fermions on the other [9]. This can be done because from (3.14), 

8 & 3 (x, y) = iu(x, 0W 3 (x, y) + 



Z (3(x,0)+ l -P(x,y) 
6 6 



& 3 (x,y) for y = 0,irR. (4.17) 



Therefore & l 3 (x, y) transforms homogeneously like the Higgs in (4.16) but without the shift 
symmetry. From the above we see we can actually provide Yukawa couplings between dif- 
ferent branes as Qii(x, 0)& t 3 (x, 7i\R)&r(x, ttR), but U(l)y invariance requires the standard 
coupling for the charge conjugates Q R and t° L to be on the same brane. We may similarly 
provide couplings for the other quark families, and also for the leptons, where we just use 
the fact that for <£p, ( |4.17| ) with y = is the usual Higgs transformation. 

Such Wilson line interactions can arise from integrating out heavy bulk fields [8, 9, 13] 
and also arise in String compactifications [9, 56]. There is no problem here and later with 
brane localised anomalies since they may be cancelled by an appropriate bulk Chern-Simons 
action [57]. However as noted by ref. [9] such Yukawa interactions produce new brane- 
localised divergences. In particular they will give quadratically divergent contributions to 
the Higgs mass 

~ AVg* 1 "^ (4.18) 

just as in the Standard Model. Of course this destroys the purpose of the model as a 
solution to the hierarchy problem. 

The reason is immediately clear because the y degree of freedom plays no role here 
for the brane localised fermions. As far as they are concerned <5 l 3 looks just like the 
Standard Model Higgs. Indeed we would also have to add brane kinetic terms and a brane- 
localised quartic potential for <3?* 3 to absorb logarithmic divergences and make the model 
r enor malizable . 1 1 

In ref. [9], the authors are interested only in the little hierarchy problem, and circum- 
vent the further difficulty (4.18) for the top by adding a new colour triplet fermion (for 



11 It is not clear to us whether the resulting framework is renormalizable in 5 dimensions, however if 
it is not, then the right approach is to consider these effects expanded around 4 dimensions, where it is 



renormalizable - see our comments below (3.11) 
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us this would be a I—1/3) so as to complete the representation Ql to a representation of 
SU(3)w- Then the summation over % in ( 4.1§| ) in fact runs over the complete represen- 
tation (i = 1,2,3) becoming just a contribution to the vacuum energy = 1). The 
divergences from the other fermions are small enough to ignore in this scenario since the 
model anyway has a low effective cutoff A and the corresponding Yukawa couplings are 
much smaller. With a Planck mass cutoff, we do not have this option, even for neutrinos 
and their very low masses, and clearly it is phenomenologically unacceptable to add hy- 
percharged but SU(2)w singlets to fill out representations of SU(3)w for the other quark 
families and the leptons. 

The appearance of divergences from such non-local operators seems to go against the 
standard lore that divergences arise from local interactions only. There is no contradiction 
however. These Wilson line interactions are local in space-time (x space). This is an 
approximation which is valid only in the limit of infinite string tension (in the case of 
String Theory compactifications) or infinite mass (in the case of massive bulk fields). In 
truth one should integrate over different locations x and x' where the ends of the Wilson 
line meet the fermions, with the displacement x — x' being weighted by a form factor of 
width 1/M, where M is the heavy scale. These Wilson line interactions are then indeed 
ultraviolet finite, quantum corrections being naturally cutoff at M. 

Since our model is supposed to be valid up to the Planck mass A however, we are forced 
to describe the dynamics that results in such smeared Wilson line interactions (unless 
M > A in which case we are back at square one) . We are left with the remaining concrete 
alternative in the literature which is indeed to generate the Yukawa interactions through 
couplings of bulk and boundary fields [10, 13]. 

We follow closely the ideas of ref. [13]. Thus we introduce a pair of bulk fermions 
^i(x,y) and ^(x,?/) (necessarily Dirac fermions since they are in 5 dimensions) with 
opposite y-parity so that a parity invariant mass term M*$>i^>2 (plus c.c.) can be built. 
Using ( |4.13 ), we choose a set of representations of SU (3) so that the Standard Model 
fermions, which exist only on the branes, can all couple to them and thus all the Standard 
Model fermions (including the neutrinos) can get masses from the effective Wilson line 
interactions that will be generated. Thus we can use a *$> a pair in the fundamental (3) 
representation so that ^ couples to Ql while \E' 3 couples to 6r. This is the minimal 
representation that will do the job. The minimal representation that yields components 
with the right hypercharge to couple to tR is the 6, where we use the 33 component to 
couple to t c L . (The 6 will also couple to Q R .) We deal with the other quark families in 
precisely the same way. Charged conjugated lepton doublets can couple first to the adjoint 
(8) representation; the charge conjugated right handed electron, muon and tau couple to 
the 333 component of a 10, and finally the right handed neutrinos can couple to the lo in 
the 8 representation ^q,s already introduced. 

As we have already discussed, something special is needed to get a low value 8/2tt ~ 
0.04 — 0.09. In ref. [58] it was shown that this is possible by adding of order 10 bulk 
fields in a mixture of fundamental and adjoint representations. They used scalars as well 
as fermions. Although scalars are problematic for us, they do not appear to be especially 
required. These authors assume that the correct value of sin 2 9w can be realised by wall 
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localised kinetic terms. As we have shown, this can indeed be achieved. Somewhat similarly, 
the authors of ref. [13] manage to reduce the minimum to 0/2tt ~ 0.096 by adding bulk 
fermions in large (rank 8) symmetric representations. They note however that this will lead 
to electroweak corrections enhanced by large group theoretical factors resulting in the scale 
at which the bulk weak coupling becomes non-perturbative being lowered (their cutoff by 
application of NDA) . Here we do not have this problem. 

For quite separate reasons we have also introduced a large number of bulk fermions. 
It would be very interesting to see if the menagerie of representations we have had to 
introduce to induce Yukawa matrices for the Standard Model fermions, also turned out to 
give values of 6 in the right range. However, we now hit a severe problem if we want to 
preserve the renormalizability of the model. The contributions of the bulk fermions to the 
bulk (3 function is 12 

3j3o = -8 T R = -4 (1 + 5 + 6 + 15) = -108 (4.19) 
R 

completely overwhelming the contribution of +11 from Am- Thus the bulk theory is no 
longer asymptotically free in four dimensions and cannot be supported by an ultraviolet 
fixed point in 4 + e dimensions. 

This constraint would appear to rule out renormalizable models based on a bulk SU(3) 
Yang-Mills theory. The obvious route to try to make further progress would be to con- 
sider larger gauge groups, thus increasing the gauge- field contribution to (5q, while also 
allowing more Standard Model fermion representations to couple to the same bulk fermion 
representations. (For example above, the 3 is shared by Ql and 6r, while the right handed 
neutrinos couple to the same representation as the lepton doublets.) This results in consid- 
ering a kind of grand unification of a quite different sort from the standard four dimensional 
cases, but with some similar properties, for example charge quantization and mass rela- 
tions. However, we also clearly need to consider more complex orbifolds, and/or deal with 
multiple Higgs vacuum expectation values arising from the Hosotani mechanism (see ref. 
[59] for such a study). 

5. Summary and conclusions 

Extra dimensional field theories have the potential to solve many of the enduring mysteries 
of theoretical particle physics. However approaches that directly address weak-scale physics 
and in particular the hierarchy problem, suffer from a severe drawback in that they are 
not renormalizable, at least as conventionally envisaged. This results in an irreducible 
uncertainty of typically ~ 1% in any predictions following from these models, and also 
implies the existence of a scale much smaller than the Planck mass (~ 100 TeV) where 
something other than field theory has to take over. 

It seems possible that a restricted class of such theories may however be made renor- 
malizable by basing their continuum limit around a non-perturbative fixed point (rather 
than the perturbative Gaussian fixed point that supports the Standard Model). 

12 The factor 8 comes from the usual 4/3, the two bulk fermions and the three families. Note that 
boundary fields and couplings make no contribution to the bulk (3 function. 
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Although it has been recognized since the 1970s [23] that non-Abelian Yang-Mills 
theory might have such extra dimensional ultraviolet fixed points, only a few studies have 
been made to search for these. From the lattice studies it seems clear that the simplest 
Wilson plaquette bare action does not allow these fixed points for SU(2) Yang-Mills in 
D = 5 or 6 dimensions. However, as we emphasised, there is no reason to expect the 
simplest action to be the correct bare action in this case. 13 The lattice study in ref. [27] 
suggests that even with some more general actions the fixed points do not exist in D = 6 
dimensions and large N. On the other hand, the exact renormalization group study by Gies 
[24] suggests that these fixed points do exist in D = 5 dimensions for SU(N) Yang-Mills 
at least for N < 5. 

Somewhat surprisingly, no-one seems to have carried the initial Wilson epsilon expan- 
sion investigation beyond the early two-loop computation of Peskin [23], so in this paper 
we do that by extending the investigation to the four loops now available. We find that 
the e expansions are very well behaved asymptotic series, for example we predict that the 
coefficients in the large N limit of the expansion for the fixed point coupling d* do not 
start to diverge until ~ 10 loops. We give values for both the fixed point coupling and 
the critical exponent u, together with estimates of the error, by following -where justified- 
the simplest methods possible. In broad terms, it seems clear that the D = 5 dimensional 
fixed points do exist, for all N > 2. The evidence for fixed points in D = 6 dimensions, is 
marginal, while we find strong evidence that the fixed points do not exist in any dimension 
D > 7. As we sought to emphasise, on the one hand e expansions above the critical dimen- 
sion such as this, are justified from studies in other models, and on the other hand there 
is considerable room for improvement on the present study by using more sophisticated 
methods. 

More generally, within an expansion in e of the D = 4 + e dimensional theory, these 
ultraviolet fixed points exist if and only if the four dimensional theory has only asymptoti- 
cally free couplings. By using a dimensional regularisation and a minimal subtraction type 
scheme, we can renormalize the theory perturbatively in the normal way in four dimensions. 
Where we need to investigate renormalization group properties in the higher dimensions, 
we can use equations such as ( |2.6D to analytically continue the results to e > 0. These 
observations open the door to constructing renormalizable extra dimensional models. 

In sec. H|, we extended this generalisation of renormalizability to include matter and 
branes. We can add fermions to the bulk providing there are not too many to destroy 
the D = 4 dimensional asymptotic freedom. It is possible to include bulk scalar fields but 
only for careful choices of couplings. The inclusion of both bulk fermions and scalars is 
only natural in supersymmetric theories (where such fixed points have been independently 
discovered [32]). Abelian gauge fields cannot however be added to the bulk. 

Several ideas in the literature such as Higgsless theories, or the inclusion of Wilson 
lines directly in the bare action to generate Yukawa terms, do not extend to renormalizable 
theories in this way, at least without further development (see sees. ^ and §). 



13 Actions with more derivatives than (2.2) would generically have problems with locality at energies of 
the overall cutoff, but this cutoff can be taken to infinity. 



When including branes, care has to be taken to work in the correct critical dimension 
for the divergences being studied. Thus in the case of D = 5 dimensions compactified 
on an orbifold S 1 /Z2, brane kinetic terms are required from renormalization in D = 5 
dimensions [47] (and not D = 4 analytically continued to 5). Furthermore, we showed that 
even if these bulk-generated divergences were eliminated in a particular model, brane-bulk 
interactions would lead to a non-renormalizable theory if the brane kinetic terms were set 
to zero. The fact that these are avoided for non-zero brane kinetic terms is a gauge theory 
analogue of the DGP effect [49], and incidentally implies that the beta functions of the 
brane-localised fields must diverge in the limit that the brane kinetic terms are turned 
off. These observations have relevance for all such extra dimensional models, not only the 
renormalizable ones being proposed here. 

We further pursue the phenomenological and theoretical constraints placed on renor- 
malizable extra dimensional models in sec. ||]. We concentrate on the weak interactions 
and gauge-Higgs unification via the Hosotani mechanism [34]: if such models can be made 
renormalizable, they become solutions to the hierarchy problem. We focus only on the 
simplest model of gauge-Higgs unification based on D = 5 dimensional SU(3) Yang-Mills 
theory compactified on S 1 /Z2. For the model to be phenomenologically acceptable, and 
perturbatively renormalizable in the manner we have described, we must have large values 
of 1/b 2 ~ 10 for the U(1)y brane kinetic terms, and a much smaller Hosotani vacuum 
angle 9 than we would find without some special mechanism. Indeed we can choose a small 
a > 0.02 (implying that the effects of the fixed point will be felt only at energies 10 times 
higher than the compactification scale, i.e. at energies ~ 20 - 50TeV) and choose natural 
values for the other brane kinetic terms, in which case we have to arrange the model to 
dynamically determine 6/2tt ~ 0.04 — 0.09. Such values would however ensure that at 
low energies the distortions of the Standard Model are sufficiently small not to come into 
conflict with precision measurements. 

The known fermions must live only on the brane and do not make a contribution to 
the bulk (3 functions (at least at one loop). The real problem arises when we consider 
how to generate effective Yukawa couplings. We seem to be forced to add bulk fermions in 
representations and number that are too large to maintain asymptotic freedom in D = 4 
dimensions, thus destroying the non-trivial ultraviolet fixed point in 4 + e dimensions. 
Note that there is no direct relation between this consequence and the fact that in four 
dimensions one cannot maintain asymptotic freedom and add sufficiently many scalars to 
spontaneously break all directions in the gauge group [44]: the group in higher dimensions 
can be much larger and broken in the first place by the orbifold boundary conditions. 

Indeed, in order to make progress, one should consider a larger group and more involved 
compactifications. As we noted in sec. [|, the Standard Model fermions can be encouraged 
to share their interactions with the bulk fermions, again reducing the problem, and leading 
to a new type of unification, the possibilities and consequences of which deserve further 
exploration. 
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